Adiabatic Approximation for weakly open systems 
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We generalize the adiabatic approximation to the case of open quantum systems, in the joint limit 
of slow change and weak open system disturbances. We show that the approximation is "physically 
reasonable" as under wide conditions it leads to a completely positive evolution, if the original master 
equation can be written on a time-dependent Lindblad form. We demonstrate the approximation 
for a non-Abelian holonomic implementation of the Hadamard gate, disturbed by a decoherence 
process. We compare the resulting approximate evolution with numerical simulations of the exact 
equation. 
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I. INTRODUCTION 

Recently, there has been a growing interest in the adia- 
batic theorem |f( in the context of quantum information, 
in particular for fault tolerant holonomic quantum com- 
putation 01, and for the design of quantum adiabatic 
algorithms 0-GI In this paper, we put forward a type of 
adiabatic approximation with focus on how an ideal (uni- 
tary) adiabatic evolution governed by a time-dependent 
Hamiltonian H is perturbed by open-system effects. We 
consider the effects of a disturbance D t /T on a desired 
ideal evolution g — ~i[H(t/T), g], as a master equation 
g = L t / T g = —i{H(t/T), g]+FD t / T (g), where T gives the 
"strength" of the disturbance, and T is the run-time. In 
the ideal case T = 0, the adiabatic approximation decou- 
ples the evolution of the instantaneous eigenspaces of H. 
fn the present approximation, the eigenspace structure 
of H(t/T) still plays a role in that it determines what 
should decouple in the adiabatic limit. It turns out that 
this in general limits the applicability to systems that are 
weakly open, i.e., to small T. The present study general- 
izes Refs. Hi, as it allows degeneracy of the Hamilto- 
nian, which is an essential feature to obtain non-Abelian 
holonomy effects in general and holonomic quantum 
computation Q in particular. 

The concept of adiabaticity in open systems has been 
addressed recently by Sarandy and Lidar in Refs. HQ for 
master equations of the above type. In their approach, 
the adiabatic approximation is characterized by a decou- 
pling in terms of instantaneous Jordan blocks of the su- 
peroperator L t / T . In other words, the decoupling is de- 
termined, in Ref. [jj, by the total superoperator L t /x, 
while this is determined by H(t/T) in the present ap- 
proach. 

The approach in Ref. Q may be difficult to use in cer- 
tain applications. One example is the analysis of holo- 
nomic quantum computation in the presence of open sys- 



tem effects. First of all it should be noted that, although 
the Jordan decomposition always exists, it could be chal- 
lenging to determine it in practice for more than a lim- 
ited class of disturbances and systems. Furthermore, for 
the approximation in Ref. || a new Jordan decomposi- 
tion has to be calculated for each choice of disturbance 
of the ideal gate. In the present approach, where the 
eigenspaces of the Hamiltonian are primary, the approxi- 
mate equation can be obtained irrespective of the form of 
the disturbance TD t /T{o)- This is due to the fact that the 
spectral decomposition of H(t/T) is in general known for 
holonomic implementations of quantum gates. For other 
applications, however, the preferred method of approx- 
imation has to be decided from the specific problem at 
hand. 

Another generalization of adiabaticity to open systems 
has been considered in > but f° r a specific type of open 
systems in the context of quantum adiabatic search. 

The structure of the paper is as follows. The approx- 
imation scheme is stated in the next section. In Sec. IIIII 
we show that the approximation can be obtained as an 
adiabatic weak open-system limit of master equations. 
Section lTvl demonstrates that the approximation leads to 
a completely positive evolution under wide conditions. In 
Sec. we apply the present approximation scheme to a 
decoherence model of a non-Abelian implementation of 
the Hadamard gate. Moreover, we compare with numeri- 
cal solutions of the exact equation. The range of applica- 
bility of the approximation is discussed in Sec. IVII The 
paper ends with the conclusions. 



II. THE APPROXIMATION 

We consider master equations of the following type 
d 



dt 



g{t) = -i[H(t/T), g(t)] + TD t/T (g(t)), (I) 
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where H{t/T) is a family of Hermitian operators, D t /T 
is a superoperator, T is the run-time of the evolution, 
and r is a strength parameter of the open system effect. 
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With the change of variables s = t/T one obtains 

±-pia) = -iT[H(s), p(s)} + TTD a (j>{8)), (2) 
ds 

where p(s) — g(sT). The superoperator D s is assumed 
to be linear. In addition to purely technical assumptions 
on D s such as sufficient smoothness with respect to s, 
we assume that the solution g(t) does not grow without 
bound with respect to some operator norm, as t grows. 
This is necessary if g(t) is to be a density operator, and is 
achieved under suitable conditions if D s is on the Lind- 
blad form. 

We assume that the dimension of each eigenspace of 
H(s) is fixed, so that we may write 

K 

H(s) = Y / E k (s)P k (s). (3) 

k=l 

Furthermore, for each s we assume Ek(s) ^ Ei(s) for all 
k, I such that k ^ I, and P k (s) are projection operators 
such that P k (s)Pi(s) = 6 M P k (s) and £ fc P k (s) = 1. 

Under conditions that are elucidated in Sec. IIIII the 
adiabatic approximation of Eq. takes the form 

p = -i[TH(s) + Q(s),p] (4) 

+rr J2 9kik'i'Ph(s)D s (p k/ {s)pP v {s))Pi{s), 

klk'l' 

where 

Q(s)=^P fe ( S )P fe (s) (5) 
ft 

is Hermitian (see Eq. 1)57^1 and g k i k 'i' are or 1 depend- 
ing on the pairwise eigenvalue differences 

A kk ,{s) = E k {s)-E k ,(s), (6) 

as is described in Sec. IIII Bl In the case of closed evolu- 
tion r = 0, we retain the standard adiabatic approxima- 
tion 

An alternative form of Eq. Q may be obtained by 
making the change of variables 

p(s)=U(s)p(s)tf( s ), (7) 

where U (s) is any sufficiently smooth family of unitary 
operators such that 

U(s)P k (s)U^(s) =P fe (0), Vfc. (8) 

In terms of p(s) , Eq. (J2J takes the form 

? = -iT[H(s), p(s)} - i[Z(s), p(s)} + TTD s (p), (9) 

where 

H(s) = U(s)H(s)UHs) =J2 E k(s)P k (0), 

k 

Z(s) = iU{s)U\s), 
D s (p) = U(s)D s (U\s)p( S )U(s))U\s), (10) 



and we have used that U(s)U^(s) = -U(s)W(s). We 
decompose the density operator as p = J2 k iP^ kl \ where 
~(ki) = p fc (o)pP,(0). We refer to p^ as the "diagonal" 
terms, while we refer to p^ kl \ with k ^ I, as the "off- 
diagonal" terms. The approximate Eq. (JIJ can be written 
as 

^ = -iTA kl (s)p( kl Hs) 
as 

-iZ k {sp kl \s)+ip^ kl \s)Z l {s) 
+TTY,gki k >iiP k (0)D s (p < - k ' l ' ) )Pi{0),{ll) 

k'V 

where Zi(s) = Pi(0)Z(s)Pi(0). The properties of g kik w 
imply that the diagonal terms p^ 11 " 1 always evolve accord- 
ing to the following equation 

£pt« = -i[Zt(s),p( u \s)} 

+TTj2Pi(0)D s (p {kk) )Pi(0). (12) 

k 

The first term on the right-hand side of Eq. 1)12(1 yields 
the non-Abelian holonomy [jj of the standard adiabatic 
approximation, while the second term introduces a cou- 
pling between the diagonal terms of the density opera- 
tor. Equation (|12|l implies that for the approximate evo- 
lution the diagonal terms always evolve independently 
of the off-diagonal terms. In the simplest case where 
gklk'V = S kk 'Sw, for k 7^ I, the off-diagonal terms evolve 
independently of each other and of the diagonal terms. 

We note that if p^ kl ^(s) are the solutions of Eq. ifTTjl. 
then W(s)p (k ^(s)U(s) = P k (s)p(s)P t (s), where p(s) is 
the solution of Eq. (@J. This follows from the fact that 
Ea. (|ll(l is equivalent to Eq. (@J, as is demonstrated in 

sec-inna 



III. THE APPROXIMATION AS AN 
ADIABATIC WEAK OPEN-SYSTEM LIMIT 

Here, we put forward one possible way to justify the 
above approximation scheme. First, we note that Eq. @ 
may be written as 

±p(s) = L^p(s) + L^p(s), (13) 

where 

= —iT[H(s), •], 

= -i[Z(s),-}+TTD s . (14) 

Since [H(s),H(s')} = 0, it follows that [L^\l^] = 0. 
This implies that Eq. 1)1 3JI can be rewritten as the follow- 
ing integral equation 

e A ^p(s) = e A(s) p(0) 

+e A ^J\- A ^L^ (e A ^p(s'))ds', 

A(s) = f S L^(s")ds", (15) 
Jo 
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where we have made the change of variables 

p(s)=e A ^p(s)- (16) 

The superoperator Lg is anti-Hermitian with respect to 
the Hilbert-Schmidt inner product (A, B) = Tr(A*B). 
Thus, exp[A(s)] is unitary, and we can rewrite Eq. 1|15|) 
as 

p(s) = p(0) 

+ J\- A ^L^{e A ^p(s'))ds'. (17) 

Note that a = exp[A(s)]c(0) is the solution of the equa- 
tion a — —iT[H(s),a], Since H(s) possesses a time- 
independent eigenbasis it follows that the corresponding 
evolution operator can be written as 

V(s) = ^exp[-zT/ fc ( S )]F fc (0), 



h(s) = / E k (s')ds' 
h 



(18) 



Thus, 



= A ( s ),r = 



a = V(s)aV\s) 

= ^e- lT{7fc(s) - /!(s)} P fe (0)aP ; (0), (19) 

kl 

for every linear operator a. We obtain 

p(s) = e- A ^p(s) 

= Y, e ' T{Ik{s) ~ I,(s)}p k(0)p{s)Pi(0). (20) 

kl 

If Eq. (|2*U|) is combined with Eq. ((T7|) the result is 

p(s) = p(0) (21) 
+ E f e iTlMk ' l ' i "' ) 



klk'l 



i JO 



xP fe (0)L^ ) (P fe ,(0)p(s , )fl'(0))P i (0)ds', 



where 



hwv{s) = 4(s) - - /fc'(s) + I v {s). (22) 
Inserting Eq. (|Tl|l into Eq. J2U yields 
p(s) = p(0) (23) 

-» E f elT4fc ' ^ ^ (°) z ( s ')^ (0), p{s')]ds' 
kk> Jo 



klk'l 



-IT ^2 I e iTI " ,h ' l ' < ' s ' ) 



, Jo 



x P k (Q)D S , (P k , (0)p( S ')Pv (0))P (O)ds', 
where we have introduced 

4z(s) =J fc («) -/,(«). (24) 



A. The diagonal terms 

The diagonal terms of Eq. 112311 read 
= P(0)p( S )P(0) 



pW(0)-i / [Z ; ( S '),P ( "V)K 



+TT V f P (0)Pv (p( fefe ) (s'))P (O)ds' 
>„ Jo 



-X rf (s). 



(25) 



Here. 

X d (s) 



= E f e lTI ' kis,) Pi(0)Z( S ')P k (0)p^( S ')d S ' 

k:k^l ^° 

- E f ^ TI ^ s y {lk) ( S ')Pk{Q)Z{ S ')Pmds' (26) 

+rr 2 r e - iT ^"( s ')p(o)p^(p( fe ' r )(s'))P(o)d s ', 

k'l'-.k'^V •'° 

where we here have used that hik'i'{s) = —I k >i>(s). 

We now show that the operator A"d(s) vanishes in suit- 
able limits of T and T. First, we cite Lemma 7.2.17 from 

0. 

Lemma 1 Suppose the function h(s) is real valued, has a 
continuous second derivative on the closed bounded inter- 
val [0, 1], and -j^h(s) ^ for all s G [0, 1]. Let the func- 
tion /(s) have a continuous derivative on [0,1]. Then, 
for sufficiently large T, there exists a constant C such 
that 



p iTh{s) 



f(s)ds < CT-\ 



(27) 



Since the integrands in Eq. Ij26|l all take the form 
exp[iTh(s)]F(s) it may be tempting to use LemmaP or 
similar results like the Riemann-Lebesgue Lemma |l3l . 
directly on these integrals. However, since F depends 
on the solution p, one should keep in mind that the so- 
lution p (and hence p) depends on T, and may contain 
fluctuations growing with T, which potentially may can- 
cel the averaging effect of the phase factors exp[iTh(s)]. 
This makes a direct use of Lemma ^ dangerous when ap- 
plied to terms containing p. In other words, we cannot 
allow the function / in Eq. (|27|l to depend on T, neither 
directly nor indirectly. It is, however, quite straightfor- 
ward to avoid this problem in the present case. 

Let {|«)}n be some fixed orthonormal basis, indepen- 
dent of s, T, and T. With respect to this basis the first 
integral in Eq. (|2(j[l can be written as 



J2W(n\ I e iT/ -( s ')(P(0)Z( s ')P fe (0)) mn pS( S ')^', 

mn " 

(28) 
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which is a sum containing integrals on the form 

e iTMs')/( s ')if(p(fcO( s '))d s '. (29) 

o 

Here, h(s) = ±[Ifc(s) — Ii(s)] for some k, I and K denotes 
a linear map from the operator to a matrix element in 
the matrix representation of it, i.e., K (•) = (n\ ■ \m) for 
some n,m. Note that the function / only depends on s, 
not on ~p or T . Similarly, the second term on the right- 
hand side of Eq. (|26|l can be written as a sum involving 
integrals of the form l|29l) . 

Now, by partial integration of Eq. (|29[) one obtains 



Rz(s) = K(p< kl \s)) / j Th ^f{s')ds' 



(30) 



JTh(s") ffJI 



f{s")ds"ds'. 



By differentiation of Eq. I|23[) . and by use of the standard 
operator norm ||er|| = sup||^,|| =1 ||cr|^)||, one finds 



ds 



^ kl \ s ') 



< A {d) + Bf'TT, 



(31) 



for some constants and B[ d ^ , where the index d sig- 
nifies the diagonal terms. From Eq. I|30|l it follows that 



\R z (s)\ < 



. (32) 



where we have used that \K(p^)\ < \\p (kl) \\ < 1, as a 
consequence of the fact that ~p(s) is a density operator. 

We note that £h(s) = E k (s) - Ei(s) for all s G [0, 1], 
which is nonzero by assumption if k ^ I. Furthermore, we 
assume that the family of Hermitian operators H(s) has 
an Hermitian continuous first derivative, which implies 
that the eigenvalues Ek{s) can be ordered, for each s, in 
such a way that they have a continuous first derivative 
(see Ref. |l4|, pp. 44-45). Thus, the second derivative of 
h(s) is continuous if H(s) has a continuous first deriva- 
tive. Moreover, the function / does not depend on T, 
and has a continuous first derivative if Z(s) and D s has. 
We may thus apply Lemma to the right-hand side of 
Eq. (|32[1 . from which it follows that there exists some 



constant C[ d ' such that 
<d) 



\Rz(s)\ < C*P(l + 4 d) ) T_1 +C { 1 a> B ( 1 a) T 



(33) 



The third integral in Eq. JUJ) may be treated in the same 
way, but including the extra factor TT, which results in 
terms Rd(s) bounded as 

\R D (s)\ < C ( 2 d) (l + A [d) )T + C ( 2 d) B {d) T 2 T. (34) 

In total, we find that the norm (or, alternatively, the 
elements in some matrix representation) of X d (s) is 
bounded as 



1 1 < 4 r 1 + B^ d) T + CfT 2 T, 



(35) 



for some constants A^ , S3 , and . 

Next, we prove that the diagonal terms of the solution 
of the exact Eq. converges to the solution of the ap- 
proximate equation of the diagonal terms, under certain 
conditions. 

The set of operators a such that J^k -Pfc(0)cr-Pfc(0) = 
<7, forms a linear subspace C of the space of all linear 
operators on TL. Define 

f d (s,a) = -i^fl(0)[Zi(»),fl(0)o-fl(0)]fl(0) 
l 

+TT Pi WD* (Pk (0)aP k (0))Pi (0). 

kl 

For er, a' € C and s g [0, 1], we have 
\\fd(s,a)- f d (s,a')\\ 

1 

+TTY / \\Pi(0)D s (Pk(0)(<y <7')^(0))fl(0)|| 

kl 

< (F^ + G^TT)\\a - a'\\, (36) 

for some constants and G^ d \ In the last inequality 
we have used that Zi(s) and D s are continuous functions 
of s and that there exist maxima of | |Z;(s)| | and 1 1 |Z? S | || = 
su P||o-||=i ll^'s( cr )l|j the latter following from [0,1] being 
a compact set. Note that the constants F^ and G^ 
can be chosen independently of T and T. Equation l|3()|) 
means that FW +G (d) rr is a Lipschitz constant for fd 
on the set [0, 1] x C. 

Suppose that /5rf( s ) is the solution of the approxi- 
mate equation for the diagonal terms, i.e., Eq. (|25l) with 
X d (s) = 0. Moreover, let 



p d (s) =5>(0)p(sW) = $> W (s), 



(37) 



where ~p(s) is the exact solution of Eq. (|23|l . We now 
intend to prove that ||7j^(s) — Pd( s )ll vanishes for all s, in 
a suitable limit. The error £, with respect to the standard 
operator norm, can be estimated as 



f(s',n(s')) - f(s',p d (s')))ds' - X d (s) 



< 





X d (s 



\\f(s',n(s'))-f(s',p d (s'))\\ds' 



< AfT- 1 +B {d) T + ci d) T 2 T 

+ (i?(<J) + G^TT) f £{s')ds'. 
Jo 



(38) 



From the above inequalities one obtains an integral in- 
equality for the error £ (s). This integral inequality can 
be shown to have the solution 



WMs)-Pd(s) 

< (A^T- 1 + B 



(39) 



(d)p + c {d) T 2 T ^ e s(F^+G^TT) 



■5 



One can conclude that a sufficient condition for conver- 
gence of the approximate and the exact solution is the 
simultaneous limits T — > oo and r — > 0, under the con- 
dition that TT is bounded. 



B. The off-diagonal terms 

The off-diagonal terms contain two types of phase 
factors, viz., exp[iTI k i(s)] and exp[iTI k i k 'i'(s)]- While 
j^^ki(s) = A k i(s) is always nonzero due to the assump- 
tion of distinct eigenvalues, the functions -j-Iklk'l'(s) — 
A k i(s) — AkH'(s) may be zero at isolated points, or more 
systematically, even if A k i(s) ^ 0. Thus, the averaging 
effect leading to the adiabatic decoupling of the exact 
equation depends upon whether the graphs of the func- 
tions Aki(s) avoid each other, cross, or coincide. We 
now study the following two physically reasonable spe- 
cial cases. 

(i) For each pair (k,l) and (k',l r ) it holds that 
A kl (s) = A kn ,(s), Vs e [0,1], or A kl (s) ^ 
A k/l ,(s), Vse[0,l]. 

(ii) For each pair (k, I) and (k 1 , 1') it holds that 
A kl (s) = A fcT (s), Vs e [0,1], or A k i(s) ^ 
A k 'i'(s) for all s, except possibly at isolated points. 
At each such point s it holds that j-(Afc/(s) — 
A^(S))^0. 

In the first case the condition says that the functions 
Afcz(s) and A k >i'(s) either coincide at all points, or never 



cross. In other words, the difference A k i(s) — A k n\s) 
is either zero or nonzero on the whole interval [0, 1]. In 
the second case we allow the above mentioned graphs to 
cross at isolated points. At those points where they do 
cross we put a restriction on how they cross, in form of 
the first derivative of A k i — A k >i>. 

As the approximation depends on the behavior of the 
functions A k i we need to keep track of those which coin- 
cide systematically. To do this in the above two special 
cases, we define 

_fl if A kl ( s ) = A k , l ,(s), Vse[0,l], 
9klk'l' -< [Q elge ^ 

Hence, g k iwv = 1 only if the two functions A k i(s) and 
A k 'i'(s) coincide systematically. Note that this definition 
holds for all combinations of k, I, k' , I', including those 
involving diagonal elements. It follows that 

9klk'V — Qk'l'kU 9klk'V — giki'k' , (41) 

and 

9kik'i = S kk ', g k iki> = 5u>, 

gkkk'v = Sk'i', 9kik'k' = S k i, (42) 

where the latter conditions hold under the assumption 
E k {s)^Ei{s),k^l. 

The equation for the off diagonal term p^> (s) can be 
written 



p(ki) {s) = p(«)( )-i / Z k ( S ')p^(s')ds' + i P^Hs'W^s' 
Jo Jo 

+TTY,9kWV f P fc (0)^(p( fe ''')(s'))^(0)^' + ^o(s), (43) 
k'V ^° 

X (s) = -i I* e lTIkk ' {s,) Pk(0)Z( S ')P k ,(Q)p( k ' l Hs')d S ' 

k'-.k'^k'' 

-i [ j TIk ' l{sl) P {kk '\s')Pk' (0)Z{a')Pi(0)da' 

k'-.k'^r 

+TT J2 / S e lT/fcifc ' i ' (s ' ) ft(0)^(ft'(0)p( S ')fl'(0))fl(0)^'. (44) 

I 



As in the proof of the approximate equation of the diago- 
nal terms we need a Lipschitz constant. Consider the lin- 
ear operators er which fulfills J2 k i-k^i -Pfe(0)crP/(0) = a. 



This operator subspace is denoted by C 1 - . Define 

fo(s, <J) = 

-i ]T Z fc ( S ')Pfc(0)aP ; (0)ds' 
ki-.k^i 

+i Y, Pk^aP^Ziis') (45) 

kl-.k^l 

+TT Y Y9kik'vPmDAPk'{o)<yPv{o))pm. 

kl:kjil k'V 
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With a similar reasoning as for the diagonal terms we 
obtain 

\\fo(s,a)-f (s,a')\\ < (F^+G^TT)\\a-<r% (46) 

for all s € [0, 1] and a, a' € £ x . Here, and G^ are 
constants. We further define 

p a (s) = ]T Pk(0)p(s)Pi(0) = E pW(s), (47) 

kl-.k^l kl-.k^l 

where ~p(s) is the solution of Eq. (|23|) . Similarly we let 

M*)= E Pa l) (s), (48) 

where pi fe ^(s) are the solutions of Eq. I|43|) . with X Q (s) = 
0. Clearly, both ~p (s) and Po(s) belong to C ± . 



1. Case (i) 

Here, the functions Aki(s) and Ak'i>(s) cither coincide 
at all points or never cross. Following the reasoning of 
the diagonal case, one finds that 



\X (s)\\<Ar'T 



(01)^-1 



B { 3 ol) T + CP'T'T (49) 



~l(ol)p2r 



and we may use the Lipschitz condition in Eq. (|46ll to 
obtain 

\\p (s)-p a (s)\\ (50) 
< (A^T- 1 + B { 3 ol) T + c { 3 ol) T 2 T)e^ F(o)+Gl ° >rT \ 

for all s e [0,1]. Thus, as for the diagonal terms, we 
find the conditions T — > 00, T — > 0, and TT bounded, for 
convergence of the approximate and exact solution. 



2. Case (it) 

In this case we allow the graphs of the functions Aui (s) 
to cross, but only at isolated points, and with a nonzero 
angle. Consider two distinct pairs (k,l) and (k',V). It 
follows that the function A^(s) — A^'i'is) has only iso- 
lated zeros. Since the interval [0, 1] is compact and since 
Aki(s) — Ak>i'(s) is continuous, there can only be a finite 
number of isolated zeros. We may partition the interval 
[0, 1] into subintervals were each subinterval has at most 
one zero of Aki(s) — Ak'i'(s) in its interior. Due to the ze- 
ros of Aki(s) — Ak'i' (s), Lemma^is no longer applicable, 
but we may instead use the stationary phase theorem, 
which we cite from 0] (Theorem 7.2.10). Note that we 
here present a weakened form of the theorem, which pre- 
cisely covers the aspects we need. 

Theorem 1 Let h(s) be analytic in a neighborhood of the 
closed bounded interval [a, b] and be real on [a, b] . Let /(s) 
have a continuous first derivative on [a, b}. Lf -j^h(s) = 



at exactly one point sq £ (a, b) and if the second deriva- 
tive of h at Sq is nonzero, then for sufficiently large T 
there exists a constant D such that 



JTh{s) 



f{s)ds < DT- 1 ' 2 . 



(51) 



We write X a (s) as a sum of integrals of the form l|29|l . 
each of which is decomposed into integrals on subinter- 
vals. In the present case, we may reason in the same way 
as in the steps from Eq. (|30fl to Eq. i|32|) . with the ex- 
ception that each integral spans only a subinterval. Note 
that we only need to use Theorem ^ on neighborhoods 
of the points where the functions Ajm cross. On the rest 
of the interval we may use Lemma ^ Thus, in order to 
use Theorem ^ the eigenvalues Ek(s) only have to be 
analytic functions of s in a neighborhood of each point 
So where Aki(so) — Ak>i'{so) — 0. Since the Hilbert space 
is finite-dimensional and since H(s) is Hermitian, this is 
the case if H(s) is analytic in a neighborhood of each Sq 
(see Ref. [ljj, PP- 33-34, or Ref. |l|, Theorem XII.3). 
We also require that Z(s) and D 8 have continuous first 
derivatives in s. 

The value of the integral 



iTh(s 



f( 3 )di 



(52) 



is 0(T 1 ) if the subinterval [a, b] does not contain a zero 
of Au(s)-A k <i'(s), and C^T- 1 / 2 ) if it does. When sum- 
ming up the contributions from the subintervals, it fol- 
lows that the value of the integral in Eq. (|32|l , for suffi- 
ciently large T can be bounded as d[ o2} T-^ 2 , for some 

constant d[ o2 \ Thus, the first two terms on the right- 
hand side of Eq. I|44|) are bounded by a finite sum of 
expressions on the form 



\Rz(s)\ < (1 + a[ o2) + b[ o2) tt)d[ o2) t 



(o2) T _l/2 



(53) 



for sufficiently large T. Similarly, for the third term 
on the right-hand side of Eq. (|44|l gives a finite sum of 
bounds of the form 

\Rd(s)\ < (1 + 4° 2) + B^T^D^TT 1 / 2 . (54) 
Thus, for sufficiently large T we obtain 

||*o(«)|| < A^T- 1 ' 2 + B { ° 2) TT 1 / 2 + C [ ° 2) T 2 T^ 2 . 

By combining this with the Lipschitz condition (|46|l . one 
obtains 

\\p (s)-K(s)\\ <{Ai o2) T-^ 2 + Bf^ 2 

+cfVT 3 /V (F<0)+G< ° )rT \(55) 

for all s e [0,1]. Since, TT 1 ' 2 = (TT)T^ 1 / 2 and 
r 2 T 3/2 _ (YTfT- 1 / 2 , it is sufficient with the simulta- 
neous conditions T — ► 00, r — > 0, and TT bounded, for 
the error to vanish. Although we obtain the same condi- 
tions as in case (i), Eq. (|55|l nevertheless indicates worse 
scaling properties of the error than Eq. (|39|l does. This 
point will be discussed further in Sec. IIII Dl 



7 



C. The approximate equations 

In Sees. IIII Al and IIII Bl we have motivated the ap- 
proximate equations for diagonal as well as off-diagonal 
terms p( kl \s). For the diagonal terms the approximate 
equation is Eq. with Xd(s) = 0. For the off- 

diagonal terms the approximate equation is Eq. 1(43(1 with 
X a (s) = 0. One may transform the integral equations 
into differential equations, followed by a change of vari- 
ables back to (s). With use of the definition of gkik'l' , 
this results in Eqs. (|12J) and 1(11(1 . for the diagonal and 
the off-diagonal terms, respectively. Note that Eq. l(TT)l 
holds, not only for the off-diagonal terms, but for the 
diagonal terms as well. This is the case since the expres- 
sion in Eq. JTTJ reduces to Eq. JT2J), due to Eq. if 
we consider the diagonal terms. 

The transformation to Eq. (0J from Eq. 1(11(1 is straight- 
forward for the dissipator. The only part which may need 
comment is the operator Q(s) in Eq. |0J. If one trans- 
forms from the variable p(s), back to the variable p{s), 
combining all the terms, one obtains 

p = -iT[H(s),p] 

+TT J2 9klk'l-Pk(s)D s (P k ,(s)pP l ,(s))P l (s) 

klk'l' 

- Pk( S )P k ( S )p - pJ2 Pk(s)P k (s). (56) 

k k 

By differentiating P|(s) = P k {s) one obtains 
Pk(s)P k (s) + Pk(s)Pk(s) = Pk{s)- If this expression 
is summed over k and is combined with the fact that 
J2k Pk( s ) = 0, the result is 

YPk(s)P k (s) = -YPk(s)P k (s). (57) 

k k 

By combining this expression with Eq. 1(56|) one obtains 
Eq. Q. 

D. Time scales 

For the diagonal terms, as well as for the off-diagonal 
terms in case (i), we have found that the error between 
the solution of the exact equation and the solution of the 
approximate equation satisfies a bound of the form 

£ < (AT' 1 +BT + CT 2 T)e F+GFT . (58) 

In view of the limiting processes considered in the pre- 
vious sections, a physical interpretation of this condition 
might be to assume that the strength parameter T de- 
pends on the run-time T. If T = a/T, with a > a 
constant independent of T, then the error would go to 
zero when T — > oo. This, however, paints our abilities to 
control open-system effects in a bit too rosy colors. In 
practice, the open-system effects are often residual un- 
controllable errors and the strength T is given by the 



situation at hand, and we have no possibility to decrease 
r as T increases. 

On the other hand, from Eq. 1)58(1 it is quite clear that 
the approximation is good if the run time T is sufficiently 
large, the characteristic time scale of the open-system ef- 
fects r _1 is sufficiently large, and the run time T is in 
the order of or smaller than r -1 . Unlike the standard 
adiabatic approximation where the error can be made 
arbitrarily small by increasing the run-time, the present 
approximation appears to be limited, since for a given 
open-system strength T, the error cannot be made arbi- 
trarily small as the run-time has to be at the same order 
or smaller than the characteristic time scale of the open- 
system effects 0| ■ One should keep in mind, though, 
that we only have obtained sufficient conditions, not nec- 
essary conditions, for the accuracy of the approximation. 
As pointed out in Sec I VII these sufficient conditions may 
in some cases be unnecessarily pessimistic. 

For the off-diagonal terms in case (ii), we similarly ob- 
tained the condition T — > oo, r — > 0, and TT bounded. 
However, this does not tell us at what rate the error 
decreases. One may compare Eqs. (|55|) and l|5%|l . As 
an example, consider those terms that solely depend on 
T, and not T. This part scales like T -1 / 2 and T" 1 in 
Eqs. lf5^|l and Eq. (jSHJ, respectively. Thus, while both 
these parts goes to zero when the run-time T increases, 
the rate is slower for case (ii) than for the diagonal terms 
and case (i). Similarly one may compare the other terms, 
in Eqs. 1(55(1 and 1(58(1 . containing combinations of T and 
r. Again one finds that the scaling of the error with in- 
creasing T and decreasing T is worse in case (ii) than for 
the diagonal terms and case (i). This suggests that the 
range of applicability of the approximation is tighter in 
case (ii). 

One may note that the constants A, B, C, G, and 
F in Eq. ((58(1 play an important role as they "set the 
scales" , in the sense that they determine what "large T" 
and "small T" means. We have avoided to give explicit 
estimates of these constants. It would be possible to per- 
form the derivations in the previous sections in such a 
way that estimates of these constants are obtained. How- 
ever, it seems a better strategy to derive such constants 
more specifically for the system and the initial conditions 
at hand. 

In essence, we have shown that there exists a region 
of large T and small T where the approximation is good, 
but we have not determined how large T and how small 
r must be. This is analogous to perturbation theory 
where one knows the approximation to be good if the 
perturbation parameter is sufficiently small, but usually 
one does not know how small "sufficiently small" is. 



IV. COMPLETE POSITIVITY 

So far we have assumed very little about the exact na- 
ture of the disturbance D s . Except that D s should be 
linear as a superoperator and be sufficiently smooth as a 



function of s, we have only required that it should lead 
to an evolution which keeps the solution p(s) bounded. 
In this section we investigate in more detail what evolu- 
tion the approximation gives rise to, and we do so for a 
restricted class of superoperators D s . 

For an important class of master equations g = 
Lg, the superoperator L can be written on the (time- 
independent) Lindblad form [l8|. If L is bounded, then 
the Lindblad form guarantees that the resulting evolu- 
tion is trace preserving and completely positive ]isl . 
To be more precise, the master equation induces a one- 
parameter family of linear maps such that p(s2) — 
A S2 _ Sl p(si), for S2 > 8\. Each is trace preserving 
and completely positive if L can be written on the Lind- 
blad form. The complete positivity guarantees that the 
evolution maps density operators to density operators, 
even if the evolution acts on one member of an entangled 
pair of systems [20j . 

If the superoperator L is time-dependent we instead 
obtain a two-parameter family of linear maps A S2jSl such 
that p(sa) = ^s 2 .siP{ s i), for s 2 > si- In the finite- 
dimensional case it can be shown that a sufficient con- 
dition for A S2)S1 to be completely positive is that the 
time-dependent superoperator L s can be written on a 
time-dependent Lindblad form, and that L s has a contin- 
uous first derivative on the interval [0, 1]. For discussions 
on the complete positivity of the dynamical map s ge ner- 
ated by time-dependent Lindbladians, see Refs. |2ll |22|. 
Note that the Lindblad form of the master equation is 
not necessary in order to obtain complete positivity, or 
positivity of the dynamical maps. More general master 
equations which are non-local in time (integro-differential 
equations) have been considered in the literature (see, 
e.g., H3)- Moreover, time- local equations p(s) = L s p(s), 
where L s is not of the lindblad form have also been con- 
sidered (see, e.g., [24|). In this section we assume that 
the disturbance D s can be written on the time-dependent 
Lindblad form. It may be possible to generalize the rea- 
soning in this section to the type of master equations 
considered in |24|. This question is, however, not treated 
here. 

Here it is shown that if the superoperator D s can be 
written on the time-dependent Lindblad form, then the 
approximate equation can also be written on the time- 
dependent Lindblad form. Thus, under suitable condi- 
tions it follows that the approximate evolution is "phys- 
ically reasonable" in the sense that it is trace preserving 
and completely positive. 

Suppose D s (p) can be written on the time-dependent 
Lindblad form 

D s (p) = -i[F(s),p}+J2Vn(s)pV^(s) (59) 

n 

n n 

where F(s) is Hermitian. To calculate the term of Eq. Q 
involving D s , we use Eq. (|59|l together with the condi- 
tions in Eq. to obtain the approximate equation 
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p = -i[TH(s) + Q(s),p] + TTR s (p), where 
R s (p) = E 9m'l>Pk(s)Ds(Pk>(s)pP v (s))Pt(s) 

Mk'V 

= -iJ2[Pk(s)F(s)P k (s),p] 

fc 

+ E 9kWi-Pk(s)V n (s)P kl (s)pP ll (s)V 1 l(s)Pi(s) 

klk'l'n 

~\j2 P k( S ) V n^)Vn( S )Pk(s)p 

kn 

-lj2p P ^ V n( S ) V n( s ) P l( s y (6°) 
In 

One may note the following 

(61) 

k 

= E ffwfcT-R'(*)V n t(s)^(fl)flfc(s)K n («)JV(*), 

klk'V 

which follows from Eq. Q42JI. By combining Eqs. (|60|) and 
(|61[1 the result is 

R s (p) = -iJ2 [Pk(s)F(s)P k (s),p] (62) 
fc 

+ E 9kwrPk(s)V n (s)P k ,(s)pPi,(s)Vj(s)Pi(s) 

klk'l'n 

-\ E 9kWvPv{s)V^{s)P l {s)P k {s)V n {s)P kl {s)p 

klk'l'n 

-\ E 9kik>i>pPl>(s)Vi(s)Pi(s)P k (s)V n (s)P k ,(s). 

klk'l'n 

Define a matrix G with elements Gkk'.w = 9kik'v ■ G is 
symmetric due to Eq. (|41|l . which implies that G is di- 

agonalizable such that G kk ',w = J2 m ^ c w c iv 4 Tnis 
can be used to show that Eq. I|62l) can be rewritten as 

R s (p) = -iJ2[Pk{s)F(s)P k (s),p] 

fc 

+ )]^MH( S )pMWt( S ) 

n m 

-\Y.Y. M( rr ] \s) M n m) {s)P 

n rn 

-^EE^ M « m)t ( s ) M ™ m) ^), (63) 

n rn 

where 

Mi m \ S ) = E V^c%)P k (s)V n (s)P k ,(s). (64) 

kk' 

Hence, we have shown that the approximate equation can 
be written on the time-dependent Lindblad form. 
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V. EXAMPLE: NON-ABELIAN HOLONOMY 

Holonomic quantum computation |2| is a recently pro- 
posed approach to quantum circuits using the idea of 
adiabatic evolution. Here, we wish to apply the present 
approximation scheme for weak open-system effects in 
holonomic single-qubit rotation gates. 

We consider a four-level system consisting of three 
ground states 0, 1, and a whose coupling to an excited 
state e is modeled by the Hamiltonian |25| 

H(s) = |e) ((OMs) + (lMs) + (o|w„(s)) + H.c. (65) 

Here, s — t/T, with T being the run-time of the process, 
and wo, u>i, and oj a are tunable, possibly complex- valued, 
coupling parameters. For each s, H possesses a doubly 
degenerate zero-energy (dark) eigensubspace spanned by 
and |x2) and two bright eigenvectors |x3) and \X4), 
the latter with energies ±w, where 

u = (66) 

This type of system is found in various implementations 
of holonomic gates, including ion-traps |26J, Josephson 
junctions [27|. semiconductor quantum dots [2^], and 
neutral atoms in cavities |29|. 

In the present investigation ui(s) is chosen to be con- 
stant. It follows that we may measure energy in units of 
u>, and thus let the vector [wo(s), ^i(s), w n (s)] to be of 
unit length. Since h = 1 it follows that we measure time, 
and especially the run-time, in units of u> . We use this 
convention in the rest of this section. 

Holonomic single-qubit rotations acting on the compu- 
tational space spanned by |0) and |1) may be obtained in 
adiabatic transport of the doubly degenerate dark states 
along paths restricted by the parametrization 

w o( s ) = sin 9(s) sin </?(s), 
LOi(s) = sin 8(s) cos </3(s), 
uj a (s) = cos0(s), (67) 

where the angles 9 and <p parametrize a 2-sphere. Ex- 
plicitly, a loop C in parameter space starting and ending 
at (u>o,u!i,Lu a ) — (0,0,1), yields the holonomic rotation 
gate 

M [ C ] = e -o(|o)(i|-|i)(o|) j (68) 

f2 being the solid angle swept by C. 

We assume that the system is influenced by an envi- 
ronment which is sensitive to whether the system is in 
the state a or not. This may be modeled by adding the 
Lindbladian 



and its concomitant strength T. 

A. Application of the approximation 

First, we notice that there is an arbitrariness in the 
choice of eigenbasis of H(s), which can be formulated as a 
choice of gauge. This arbitrariness in the choice of gauge 
is related to the arbitrariness in the choice of U(s) in 
Eq. Let {\xk(s))}k be an instantaneous orthonormal 
eigenbasis of H(s). Given such a basis one may construct 
a family U(s) by 

U(s) = U J2\Xk(0))(Xk(s)\, (70) 
fc 

where Uq is a fixed unitary operator such that 
[Uq, -Pn(O)] = for all n. Every family U(s) constructed 
via Eq. (|70|l is unitary and satisfies Eq. JSJ. Moreover, ev- 
ery family U (s) that satisfies Eq. © can be reached via 
Eq. (|70(l for some choice of instantaneous orthonormal 
eigenbasis {|xfc(s))}fc and Uq. As the present approxi- 
mation is independent of the choice of allowed C/(s), it 
follows that the approximation is also independent of the 
choice of gauge. 

Here we briefly describe a procedure to put the approx- 
imate master equation into matrix form. As the present 
application only regards the computational subspace we 
disregard the off-diagonal terms of the approximate so- 
lution, and only consider Eq. (|12|l . In the present case 
there are three diagonal terms, corresponding to the dark 
subspace and the two bright states. In order to write 
Eq. H12|) on matrix form we first choose an instantaneous 
orthonormal eigenbasis {\xk(s))}k of -ff(s), from which 
one can construct U(s) via Eq. (|70|l . with Uq = 1. Define 



P a = P12 P21 P22 Ph Ph)\ (71) 

where p^(s) = {Xk(0)\p a (s)\xi(0)), and where p a (s) is 
the solution of Eq. 112fl . Note that we here use the initial 
eigenbasis {|xfe(0))}fc. This is related to the fact that 
Eq. (|12fl is written in the "rotated frame", as described 
by Eqs. J7J) and JSJ. When inserting U(s) into Eq. ifT^I 
the result can be written as 



p a = M a {s)p a . (72) 



V=\a)(a\ 



(69) 



If the instantaneous eigenbasis is chosen to be 
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|xi(s)) = cos<p(s)|0) -sin<p(s)|l), 

IX2(s)) = sinc/?(s) cos#(s)|0) + cos <p(s) cos0(s)|l) — sin#(s)|a), 

|x 3 (s)) = —=(ma(p(s)BWi9(s)\0) + cos tp(s) sin 9 (s)\l) + cos9(s)\a) + \e) 
v2 v - 

|X4(s)) = —/=( sinip(s) sin#(s)|0) + cosip(s) sin0(s)|l) + cos9(s)\a) — \e) 
v2 ^ 

then, with the Hamiltonian in Eq. (|65fl and the Lindbladian in Eq. (|69[) . one obtains 

/ cos 9(s) cos 9{s) 

E£ sin 2 9(s) -^cos0(s) 

M a = ^cos9(s) -££sin 2 0(s) -^cos9(s) 

km 



cos 0(a) 

Q 







%cos8(s) % cos 9(s) 



§/(*) 



o 




§/(*) 

rT cos 4 0( S ) 



/(s ) ^cos 4 0( s ) -\g(s) J 



(73) 



(74) 



Here, f(s) = TT sin 2 6»(s) cos 2 6»(s) and gr(s) 
17 [1 + sin 2 (9(s)] cos 2 6»(s). 

From the above analysis it follows that the so- 
lutions of Eq. (|12f> can be written p^ nn '(s) = 
Efci Pm( s )Ix/c(0))(xz(0)|, where the sum over fc,Z spans 
the appropriate elements for each n. Since these oper- 
ators are written in the "rotated frame" , it is appropri- 
ate to invert this transformation to more easily analyze 
the gate operation. By inverting the transformation in 
Eq. (0 one finds that these operators can be written 

U\s)ft nn Hs)U(s) = 52p a kl (s)\ X k(s)){xi(s)\ 
kl 

= P n (s)p(s)P n (s), (75) 

where p(s) is the solution of Eq. (JIJ. 

There are some subtleties associated with the choice of 
basis and the usual difficulty with spherical coordinates, 
viz., that (p is not defined at the north and south pole 
of parameter space. In fact, with the choice of basis in 
Eq. (|73p . the gauge potential Z(s) has singularities at 
both poles. Nevertheless, if we avoid loops around the 
poles and take appropriate limits if we wish to approach 
the poles, then this gauge is unproblematic. Another pos- 
sibility is to rotate the dark instantaneous eigenvectors 
as 

2 

\xk)^\x'k) = El^>fe|wix fc ), fe = l>2, 

i=i 

W = e p(lxi><3C2|-|x2>(xd)_ (76) 

With this basis one obtains a gauge where the vector 
potential is well defined except at the south pole of pa- 
rameter space. The rotated dark states does, however, 
give a system of differential equations too extensive to 
be presented explicitly here. 



We restrict the paramctrization to 

<p(s) =as + b, 9(s)=cs + d, a,b,c,deR, (77) 

and the paths in parameter space to half "orange slices" 

(p = O,0 = CI,t = O) -> (0,7r/2, Ti) 

-» (Sip, tt/2,T 2 + Ti) 
-> (^O.Tg + Ta + Ti) 

-» (o,o,T4 + r 3 + T 2 + r 1 ), 

(78) 

where !Z\, . . . , T4 are the run-times for the path segments, 
see Fig. ^ Note that the fourth path originates from 
the deformation of a well defined square to the orange 
slice on the parameter sphere. In the limit — > the 
fourth path is reduced to a single point at the north pole. 
This implies that T4 can be set to zero without loss of 
adiabaticity. Note that in the present decoherence model 
a nonzero T4 only affects the two bright states. 

For the initial state vector |^) = cos(x/2)|0) + 
e~ ly sin(x/2)|l), the output state of the approximation 
projected onto the computational space may be written 
as 

p a out = u[C}p'uHC], 
P' = ( l i + i C08X n J 1 w ), (79) 

where 

h = ie-^ T3+2T2+Tl V^smx, (80) 
h = l+ 2 -e-^ T ^\ (81) 

and u[C] is the holonomy Eq. (|68|) in the {|0), |1)} basis. 
Thus, the output is determined by the holonomy trans- 
formation of the (^-independent p' . It is worth to point 
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FIG. 1: Path in parameter space starting and ending at P. 
The duration of the path segments P— >Q, Q— >R, and R^P 
is Ti, T2, and T3, respectively, as described in Eq. 1781 . 6<p is 
the opening angle in the equatorial plane. Note that for this 
path the enclosed solid angle equals Sip. 



out that this feature is due to the particular choice of 
parametrization and the chosen loop in parameter space, 
and not some intrinsic property of the decoherence. Fur- 
thermore, in addition to destroying the superpositions 
between the computational states the decoherence also 
gives an intensity loss from the computational space. 
This intensity loss arises as the states corresponding to 
X2(s), X3( s )j an( i X4( s ) decoheres into mixed states, as 
all three of them contain the a state. 

As a final observation we note that Eq. H72[) can be 
obtained more or less directly from Eq. J2J. We may rep- 
resent Eq. using an instantaneous orthonormal eigen- 
basis {\xk{s))}k as 

P = M(s)p, (82) 
where M(s) is a 16 x 16 matrix, and where 

P=(/>n P12 ■■■ p\i P21 ■■■ Pa) 1 (83) 

with pki(s) = (Xk(s)\p{s)\xi(s)) , k,l = 1,...,4. Due 
to Eq. J7J) it follows that Eq. is also obtained if 

we instead represent Eq. © using the {|xfc(0))}fc ba- 
sis, where we again assume that U(s) is constructed via 
Eq. lfrtJ|) with Uq = 1. Equation ((TTJl is obtained by re- 
moving couplings from Eq. 10. For the chosen basis, 
this corresponds to a removal of off-diagonal elements in 
M(s), such that the new approximate matrix can be ar- 
ranged in a block diagonal form. Each of these diagonal 
blocks corresponds to a collection of coupled terms. One 
block corresponds to the diagonal terms, and there is one 
block for each collection of off-diagonal terms that cou- 
ples among themselves, as determined by gkik'v- If one 
is interested in the evolution of a particular collection 
of coupled terms, then the approximate equation is ob- 
tained if one removes those rows and columns from M(s) 



that correspond to terms not included in the collection. 
In the present example, the matrix M a (s) in Eq. H74|) 
is obtained if we use the basis in Eq. I|73[) to represent 
the exact master equation, and remove those rows and 
columns from M(s) that correspond to the off-diagonal 
terms. 



B. Numerical analysis 

We compare the approximate solution in Eq. (|79|l with 
a numerical solution of Eq. 18211 in the Hadamard case 
n = 7r/4 by putting Sip = ir/4. For the calculation we 
have used the gauge where the vector potential is well 
defined at the north pole. We further put T4 = 0. We 
distribute the run-time T proportionally to the length of 
the three circle segments, i.e., I\ = T 3 = 2T/5, T 2 = 
T/5. 

In the numerical treatment of the evolution we de- 
compose the interval [0, T] into subintervals with step 
size At, on which M(t) is taken to be constant. The 
resulting approximate evolution is on the form pK — 
[nf =0 exp(AiM(£ fc ))]p - The step size At = 0.01 is 
used. The step size At = 0.005 has been tested, without 
any significant change of the result. 

For quantum gates, the relevant error is that at the 
end-point. This error may contain a contribution in form 
of an intensity loss out of the computational subspace. To 
detect this intensity loss, we use the quantities 



I(T) = l-Tv(Pg(T)), 
I a (T) = 1 -Tr (Pg a (T)) 



(84) 



with P the projector onto the computational subspace 
spanned by |0) and |1). The error within this subspace 
is analyzed in terms of the fidelity [3(], 

£>( gnor m (r), g g orm (T)) 

= Tr^ A / ftl orm(r)es orm (r)v / ^5?y, (85) 



where 



Pnorm (^J — pnorm (1) 
Pnorm 

(T) = 

Pnorm (i) 



PQ(T)P 
Tr(Pg(T)) ' 

Pg a (T)P 
Tr^P^T))"' 



(86) 



are the normalized outputs of the exact and the approx- 
imate evolution, respectively. This normalization may 
correspond to a post-selection procedure. 

In Fig. we show (l|e(T)|l) and (l\g a (T)\l), for 
r = 0,0.01,0.1. We have chosen the initial state vec- 
tor |*) = cos(x/2)|0) + e - *»sin(a;/2)|l) with x = tt/5 
and y — 3n/ '4. The corresponding normalized fidelity D 
is shown in Fig. Inland the intensity losses I(T) and I a (T) 
are shown in Fig. 0] These simulations indicate that for 
this model system the error seems to decrease with in- 
creasing run-times T at a rate more or less equal to the 
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ordinary adiabatic approximation in the closed case, in- 
dependent of the strength T of the decoherence process. 
We have confirmed this finding for other input states. 




50 100 150 200 



FIG. 2: The solid lines show the value of the matrix ele- 
ment (l|p(T)[l) of the output density operator of the exact 
equation, as a function of T. The run-time T is measured 
in units of ui^ 1 , defined in Eq. IRjfcil . The dashed lines show 
the corresponding value (l|f} a (T)|l). Counted from the top 
and down, the solid-dashed line pairs correspond to T = 0, 
r = O.Of, and F = 0.1, respectively. The initial state is pure, 
with polar angle x = n/5 and azimuthal angle y = 3n/4 on 
the Bloch sphere. The horizontal dashed line corresponds to 
the ordinary adiabatic approximation for the closed evolution 
case. One may note that the rate at which the exact solution 
approaches the approximate solution appears to be rather in- 
dependent of the strength F of the decoherence. 



VI. RANGE OF APPLICABILITY 

The analysis in Sec. lIIII suggests that, for a given T, the 
error bound in Eq. H58J1 has a minimum for some value of 
the run-time T. It is quite straightforward to obtain an 
example of a system where this appears to be the case. 
One may consider a time-dependent Hamiltonian of the 
form 

H(a) = e- lsZ H Q e lsZ , (87) 

where Hq and Z are fixed Hermitian operators. The spec- 
trum of this Hamiltonian is fixed, but the eigenbasis ro- 
tates. One may consider the master equation 

p = -iT[H(s),p]-TT[A,[A,p]}, (88) 

where A is a fixed Hermitian operator. The double com- 
mutator in the above equation causes decoherence with 
respect to the eigenbasis of A. We have chosen a four di- 
mensional Hilbert space and have generated Hq, Z, and 
A, as well as the pure initial state, randomly. Figure 
shows the maximum error in the Hilbert-Schmidt norm 
max sg [o.i] \\p{s) — p a (s)\\ for various choices of T and T. 




0.7 V 



20 40 60 80 100 

T 

FIG. 3: These graphs highlight another aspect of the same 
series of calculations as in Fig.H They show the error between 
the exact evolution and the approximate evolution, in form 
of the normalized fidelity D defined in Eqs. I85H and 1861 . as 
a function of T. The run-time T is given in units of uj' 1 , 
defined in Eq. I66H . Note the different range of the run-time 
compared to the other plots. Here the dotted line corresponds 
to r = 0, the dashed line to F = 0.01 and the solid line to 
r = 0.1. Note that the dotted and the dashed lines almost 
coincide. These graphs indicate that the distance between the 
approximate and exact evolution, within the computational 
subspace, decreases with the run-time T at a rate more or less 
independent of V. 



As seen in Fig. we indeed seem to have the expected 
behavior of the approximation. In the ideal case, T = 0, 
the error appears to go to zero as T increases, while for 
non- vanishing T there seems to be a minimum error. 

However, the error does not always seem to behave 
in this manner. In the example presented in Sec. IV Bl 
there is no trace of this minimum error. Rather the error 
seems to vanish for large T for any value of T. In other 
words, the error bounds derived in Sec. IIIII appears to 
be unnecessarily pessimistic in some cases. We here put 
forward some reasons why this may be the case. 

One aspect is the question of which error to consider. 
In Sec. lIIII we considered the maximum deviation between 
the exact and approximate solution during the whole evo- 
lution, while in Sec. IV Bl the relevant error was taken at 
the end of the evolution. In some cases the maximum de- 
viation need not occur at the end of the evolution, which 
may cause the "end point error" to be smaller than the 
maximum deviation. One may also note that Sec. IV Bl 
focused on one single diagonal term of the total density 
operator and that the error for this part may be smaller 
than the total error. 

Another reason for the approximation to be accurate 
under wider conditions is if the dissipator D s is such that 
it does not couple off-diagonal terms to diagonal terms, 
or off-diagonal terms to other off-diagonal terms. Un- 
der such conditions the dissipator is unaffected by the 
approximation and it seems reasonable that the approxi- 
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FIG. 4: The intensity losses, as defined in Eq. (1841 . out of the 
computational subspace for the same series of calculations as 
in Figs. |21 and |3] The solid lines show the intensity losses I 
of the exact evolution, and the dashed the intensity losses I a 
of the approximate evolution, as a function of T. The run- 
time T is given in units of defined in Eq. (16611 . Here 
the lowermost pair of curves correspond to V = 0. Note that 
for r = the loss is identically zero for the approximate 
evolution. The uppermost pair of curves corresponds to V = 
0.1, and the pair in the middle corresponds to T = 0.01. 



mation should have a wider range of applicability. How- 
ever, this cannot be the sole reason, as is indicated by the 
results in Sec. IV Bl since the dissipator used (i.e., deco- 
herence with pointer state a) does belong to the class of 
dissipators that do couple the diagonal and off-diagonal 
terms. Suppose, however, that the evolution is such that 
the magnitude of the off-diagonal terms tends to decrease 
with increasing run-times. For example, this may oc- 
cur if a decoherence or relaxation process acts suitably 
in relation to the instantaneous eigenspaces. Consider 
the diagonal terms: even if there would be a coupling 
to the off-diagonal terms, the importance of this cou- 
pling naturally diminishes if the off-diagonal terms tends 
to decrease in magnitude. If the open system process 
is such that it tends to suppress the off-diagonal terms, 
it thus seems reasonable to expect that the approximate 
equation should be accurate at large run-times. This re- 
duction of off-diagonal terms reasonably should be more 
relevant for the end point error than for the maximum de- 
viation. Another reason for the end-point error to vanish 
is if the approximate and exact equations have a common 
asymptotic state. 

There is clearly room for further investigations of when 
and why the present approximation is accurate beyond 
the joint limit of slow change and weak open system ef- 
fects. 
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FIG. 5: The maximum error in the Hilbert-Schmidt norm 
max s6 r 01 i ||p(s) — p a (s)|| between the solution p(s) ofEq. 1881 
and the solution p a (s) of the approximate equation as a func- 
tion of the run-time T, the latter measured in arbitrary units. 
The plots are generated for one random instance of Ho, Z, 
A, and initial state, for a four dimensional Hilbert space. 
Each curve corresponds to a value of F, and shows the maxi- 
mum error as a function of T. To the right of the figure the 
curves correspond to, counted from the bottom and up, to 
T = 0,0.002,0.004,0.006,0.008,0.01. As seen, the error for 
the adiabatic approximation in the closed case (F = 0) seems 
to tend to zero as T increases, while the other cases appear 
to have a minimum error for a certain T. 



VII. CONCLUSIONS 

We present an adiabatic approximation scheme for 
weakly open systems. Contrary to the adiabatic approx- 
imation for closed systems, the presence of open system 
effects introduces a coupling between the instantaneous 
eigenspaces of the time-dependent, possibly degenerate, 
Hamiltonian. We show that the present approximation 
can be obtained as a slow-change weak open-system limit, 
in the sense that the time scale inversely proportional to 
the strength of the open system effect puts an upper limit 
on the run-time. In the ideal case of closed systems this 
limiting time scale becomes infinite, and the ordinary 
adiabatic approximation is retained. 

We demonstrate the approximation scheme for a non- 
Abelian holonomic implementation of a Hadamard gate, 
exposed to a decoherence process. We compare the ap- 
proximation with numerically obtained solutions of the 
exact master equation. These calculations indicate that 
the error between the approximate and the exact evolu- 
tion decreases with increasing run-time at a rate more or 
less independent of the strength of the decoherence pro- 
cess. This result suggests that the approximation scheme 
may have a wider range of applicability than the weak 
open-system limit. 
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